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Abstract: We investigate on the nonlinear dynamics of Airy beams 

in a regime where nonlinear losses due to multi-photon absorption are 
significant. We identify the nonlinear Airy beam (NAB) that preserves the 
amplitude of the inward Hankel component as an attractor of the dynamics. 
This attractor governs also the dynamics of finite-power (apodized) Airy 
beams, irrespective of the location of the entrance plane in the medium 
with respect to the Airy waist plane. A soft (linear) input long before the 
waist, however, strongly speeds up NAB formation and its persistence as a 
quasi-stationary beam in comparison to an abrupt input at the Airy waist 
plane, and promotes the formation of a new type of highly dissipative, fully 
nonlinear Airy beam not described so far. 
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1. Introduction 

The discovery by Berry and Balazs JT] in Quantum Mechanics of the existence of non¬ 
spreading, self-accelerating Airy wave packets hibernated for decades in the literature until 
its experimental realization in the form of an Airy beam of light (2j[3l. It was indeed in Optics 
where the land was prepared thanks to the re-discovery in the eighties of non-diffracting Bessel 
beams 0. Since then huge amount of work has been made on self-accelerating Airy beams, 
including their properties 0161, propagation dynamics in different nonlinear media 1714 1 2L 
nonparaxial |fl3l , vectorial |jT4| and fully nonlinear fl5l generalizations, and applications as 
diverse as particle clearing Cl or creation of curved plasma channels E). It seems that this 
amount of work was necessary to close the circle with the recent creation of genuine quantum 
Airy wave packets of free electrons US). 

In this paper we are interested in the nonlinear propagation of one-dimensional Airy beams 
in a regime where nonlinear losses (NLLs) due to multiphoton absorption are significant. In this 
regime, Lotti el al El have demonstrated theoretically the existence of Airy-like beams that 
can propagate stationarily, overcoming the attenuating effects of NLLs. These nonlinear Airy 
beams (NABs) cannot have arbitrary intensities, but their peak intensity is always lower than a 
maximum value determined by the optical properties of medium. The mechanism of stationarity 
of NABs is similar to that described earlier for nonlinear Bessel beams in the same regime 
HI 911201 . Asymptotically, a NAB behaves as a linear Airy beam except that the amplitudes, say 
|a; n | and |a ou t|, of its inward and outward Hankel components are not equal. The excess of the 
inward Hankel component creates an inward power flux that refills the power absorbed mainly 
in the more prominent, central intensity maxima. NABs are characterized by a loss of contrast 
and compression (in self-focusing Kerr media) of the intensity ripples in their intensity profile 
compared to those of the linear Airy beam El- The numerical simulations and experiments 
in El suggest that a linear Airy beam introduced in a nonlinear medium with NLLs transforms 



into a stationary NAB because the loss of contrast and compression characteristic of NABs were 
indeed observed. Also, more detailed numerical simulations in ED evidence the trend of the 
input linear Airy beam towards a stationary NAB. In case that a stationary NAB is (or tends to 
be) formed, a relevant question is which particular one, among all supported by the medium, 
is selected as the final stationary state for a given input Airy beam, so that its intensity can be 
predicted in practice. Further, finite-power Airy beams, e. g., exponential-Airy beams, have a 
transversal intensity profile changing from Gaussian at the far field to exponential-Airy at the 
“waist” El- Another relevant question is then whether these different entrance conditions of the 
input Airy beam may lead to NAB formation or not, and in the affirmative case, whether the 
formed NABs are different or not. 

Our numerical simulations with different values of the intensity, apodization parameter and 
entrance plane of the input Airy beam shows in fact seemingly disparate behaviors. We show 
however that all them can be understood from a simple conservation principle. When an Airy 
beam («i n = a ou t = a o), apodized or not, at its waist or not, is launched into the nonlinear 
medium, the NAB that preserves the amplitude of the asymptotic inward Hankel component, 

1. e., that with c/ m = ao, tends to be spontaneously formed. A similar conservation principle 
has been given recently for Bessel beams EOl . The nonlinear dynamics of Airy beams is then 
determined by the asymptotic behavior of NABs, which is analyzed in detail in Sec. [2] 

With ideal Airy beams, the stationary attracting NAB is fully formed, as seen in Sec. [3] With 
finite-power Airy beams, the situation is more complex, through the trend towards the same 
attracting NAB continues to dominate the dynamics, as seen in Sec. [4] As two representative 
cases with seemingly opposite behavior, we consider exponentially apodized Airy beams en¬ 
tering the medium at their waist, where they are fully formed and more intense (nonlinear or 
“abrupt” input), and entering the medium long before the waist, where they are widespread and 
Gaussian-like, so that nonlinearities are initially negligible (linear or “soft” input). It is note¬ 
worthy that the same NAB tends to be formed under such a disparate input conditions. With 
abrupt input, however, we only observe the first steps, mimicking the ideal case, towards the 
formation of the attracting NAB, but the beam decays well before it can be reached. Actually, 
no vestige of stationary behavior can be identified. With soft input of the same Airy beam, 
the attracting NAB is found to be much more closely reached. Its transversal profile is clearly 
observable and is seen to behave as a quasi-stationary light beam. 

Also, we identify a propagation regime in the limit of very high power of the input linear Airy 
beam for which NABs with asymptotic Hankel tail are no longer formed. What is observed is 
instead a finite-power version of the limiting NAB (LNAB) supported by the medium. This 
LNAB does not have a linear Hankel tail, but is nonlinear everywhere, strongly dissipative 
(infinitely dissipative, ideally) and more weakly localized than NABs. This attracting LNAB in 
this limit of high powers is solely determined by the optical properties of the medium. 

2. Asymptotic behavior of ideal nonlinear Airy beams 

Our starting point is the nonlinear Schrodinger equation (NLSE) 

d-A = — <3 2 A + |A| 2 A - ^|A| 2M “ 2 A, (1) 

2k x n 1 1 2 1 1 ’ W 

for the complex envelope A (x,z) of the light beam E = Aexp[— i(a>t — fc)] of angular frequency 
CO and linear propagation constant k = ( co/c)n . In the above equations n is the linear refractive 
index, c the speed of light in vacuum, «2 is the nonlinear refractive index, and /j' M) >0 is 
the M-photon absorption coefficient responsible for NLLs. For a more comprehensive analysis 
we introduce the dimensionless coordinates t, =x/xo, C = z/(fc*o)> where xo is an arbitrary 
transversal scale, the dimensionless envelope A = /2) 1 ^ 2 m_ 2 'A, and describe the evo- 



lution of the envelope in the uniformly accelerated reference frame defined by w = C (£/2) 2 , 
v = With these changes Eq. ([]} becomes 

d r A = l -dlA + V -d u A + ia\A\ 2 A - |A| 2M ~ 2 A, (2) 

where a = {k 2 x^H2/n)(2/kx^^ M ^)^/ M ~ l \ In the absence of the nonlinear terms, Eq. 0 
admits the linear Airy beam solution |T| A(u,v) At(«)exp[/<jc(w,v)], where 0 l(k,v) = 
uv/2 + v 3 /24, whose intensity profile is invariant in the accelerated reference system. 
With the nonlinear terms, we search, as in CD, for stationary solutions of the form A = 
a(u) exp[; i//(m)] exp[i(j>r(u.v)], where a(u) is the real amplitude and is a nonlinear phase. 
The NLSE in Eq. © then yields 

a" = ua + { fit ') 2 a — 2aa 3 , 

yr" = _2^-2a m - 2 , ( 3 ) 

a 

where primes stand for d/du. The second equation is readily seen to be equivalent to 

yr'a 2 = 2 / a 2M du = N u , (4) 

JU 

which is the refilling condition for stationarity with NLLs, stating that the NLLs N u in [k,°°) 
must be compensated by a positive flux y/'a 2 at u. NABs are solutions of these equations be¬ 
having as a(u) —> aoAi(u), yr{u) —> 0, i. e., as the linear Airy beam, for 11 —» +°°, and verifying 
the localization condition a(u) —> 0 for u —>• —°°. An asymptotic analysis reveals several possi¬ 
bilities for the decay for 11 —> — °°, depending on the admissible values of the total NLLs N- 00 
that the beam experience per unit propagation length. 


2.1. Nonlinear Airy beams with finite nonlinear losses 


As in CD, we first assume that N -00 < Replacing N u with N in Eq. © at large enough 
negative u, i. e., y/'a 2 ^ 2 N-«,, and with the change a = b(p)/(—u )*/ 4 , where p = (2/3 )(—m) 3 / 2 , 
we derive the Newton equation d 2 b/dp 2 = -b + N 2 m /l:fi (after neglecting subleading terms in 
p), which is the equation of motion of a particle subjected to a conservative force with potential 
V = b 2 /2+ \/2N 2 i 00 /b 2 , and therefore with a conserved “energy” E = (1/2 )(db/dp) 2 + V. 
The general solution of the Newton equation leads to the possible asymptotic behavior for the 
intensity 
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where C = (1 — N 2 , X JE 2 ) 1/2 and N < E. This is the asymptotic behavior of the NABs re¬ 
ported in CD, decaying as the linear Airy beam for u —» —but with possibly shifted oscilla¬ 
tions by a phase (j> and reduced contrast C < 1. These NABs are therefore asymptotically linear 
waves, whose NLLs are indeed finite and concentrated in the nonlinear core close to u = 0. 
Considering the asymptotic expressions of Hankel functions of order 1/3 for large negative 
argument l22l . and as pointed out in CD, the above asymptotic behavior of the intensity a 2 and 
the flux yt’a 2 22 AC x is reproduced by the expression l23l 
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provided that |a ou t| 2 = 2n(E |fli n | 2 = 2n(E +AC„,), and arg(a 0 ut) — arg(aj n ) = <j). 

The linear Airy beam A = aoAi(it)exp[i$j,(K,v)] is recovered with a m = fl G ut = «o- F° r NABs 
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Fig. 1. (a) Gray curves: Intensities profiles of NABs with ciq = 1.500, oq = 4.000 and 
«0 = 4.325 in a medium with M = 5 and Kerr nonlinearity a = 1. Black curve: Intensity 
profile of the LNAB with a o = flo.max — 4.343. (b) The same but in double logarithmic 
scale. The curves decay as (— u)~ x ! 2 for the NABs and as (—for the LNAB 
for u —» —oo. (c) Nonlinear losses ALoo and oscillation contrast C = (1 — N^/E 2 ) 1 / 2 of 
NABs as functions of «o- (d) Amplitudes |ai n | and |a out | of the inward and outward Hankel 
components of NABs as functions of oq. In (b) and (c), M = 5, and a = 0 (green curves), 
and a = 1 (black curves). 


1i n |“ > |a out | 2 , indeed |a; n | 2 — |a out | 2 = 4;rfV_oo > 0. Using the relations between Hankel and 
Airy functions of the first and second type ll22l . Eq. <[6]) can be interestingly written in more 
compact form as A m pAi(u) + iqBi(u), where p = (aj n + a ou t)/2 and q = (a; n — a Q u t)/2. 

Numerical integration of Eqs. (OJ with initial condition a(uo) = a^Aifno) and (// = 0 at some 
large positive uq, shows that NABs with the above asymptotic behavior exist only for ao < 
ao.max, where ao.max is a number of the order of unity that depends on M and the strength of 
Kerr nonlinearity a (e. g., for M = 5, ao.max — 1.72 for a = 0 and ao. max — 4.34 for a = 1). 
The three gray curves in Fig. [lja) represent the intensity profiles of NABs with increasing ao 
in a medium with M = 5 and a = 1. As aq increases, the peak intensity increases at first, but 
saturates as ciq approaches ao.max. showing that NABs can exist only up to a maximum value 
of the peak intensity. The same NABs are represented in Fig. |TJb) in double logarithmic scale 
in order to appreciate their asymptotic decay as l/(— m) 1 / 2 . Figure [TJc) shows that the total 
NLLs A-,*, of NABs are finite for any ao, but grow without limit as a o —>• ao, m ax. and that 
the contrast of the oscillations decreases down to zero in this limit. The values of /V rx , and C 
allow us to find the values of the amplitudes |aj n | and |a out | of the inward and outward Hankel 
components, which are represented as functions of ao in Fig. [TJ d). These curves turn out to be 
essential for understanding the nonlinear dynamics of Airy beams, either ideal or their finite- 
energy realizations. In absence of Kerr nonlinearity (a = 0), a m and a out turn out numerically to 
be real, while with dispersive nonlinearities as Kerr nonlinearity arg(a out ) = —arg(a; n ). Similar 
asymptotic analysis has been made recently for vortex Bessel beams f20l . A particular NAB 
and its asymptotic form in Eq. © can be seen in Fig. [3 a). 

2.2. Nonlinear Airy beams with infinite nonlinear losses 

We point out that finiteness of is a reasonable but not necessary assumption. Localized 
wave objects with infinite power and infinite nonlinear losses have been previously described 
and have been shown to play a role in actual physical processes, as in the self-focusing of 
Gaussian wave packets halted by nonlinear losses ll24l . 

Let us assume a more general decay a ~ b(—u)~ a for u —» — °°, where b and a are positive 
constants, and that a is small enough so that N„ —> °° for u —> —oo. Integrating Eq. ©, we get 
t)/ ^ 2b 1M ~ 2 {— m) 1_2ci * m_1 )/(l — 2Mo), where we have neglected a finite constant compared 
to the assumed infinitely large NLLs. Introducing this expression into the first of Eqs. © we 
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Fig. 2. In a medium with M = 5 and Kerr nonlinearity a = 1. in double logarithmic scale, 
(a) intensity profile of the NAB with ao = 4.0 (solid curve), and its asymptotic form in Eq. 
© with |a; n | = 3.176 and |a out | = 0.523 (dashed curve); (b) intensity profile of the LNAB 
with ao.max = 4.343 (solid curve) and its asymptotic form a 2 ~ b 2 (—u)~ 20 with b and a 
in Eq. GQ (dashed curve). 
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for u —>■ —The second and third terms in the left hand side cancel each other asymptotically, 
and the other two terms are negligibly small, if 
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With this value of <7, N~«, tends indeed to infinity. 

From a purely numerical approach, it turns out that the transversal profiles of NABs with 
a o approaching ao.max [the upper gray curve in Figs. |TJa) and |TJb)] tend to match a certain 
transversal profile [the black curve in Figs. [Ha) andQJb)], whose behavior is described by Eqs. 
©. The matching region expands toward negative values of u as ao increases up to the limiting 
value ao.max^ and fills the whole space down to u = — °° in the limit a o = ao.max- The asymptotic 
Airy-like behavior of NABs, and in particular its structure as two interfering Hankel beams, thus 
disappears in this limit, and is replaced by the asymptotic behavior of this limiting nonlinear 
Airy beam (LNAB) with ao = ao. max . Contrary to NABs, NLLs are not concentrated in the 
nonlinear core, but are significant in the whole transversal profile. We note that this LNAB has 
no free parameters, but is solely determined by the properties of the medium (M and a in our 
dimensionless variables). Figure|2jb) shows the LNAB in a medium with M = 5 and a = 1 and 
its asymptotic behavior determined by Eq. ©. Despite the ideal behavior of the LNAB, it will 
seen to act as an attractor in the nonlinear propagation of very powerful, real Airy beams. In 
another context, fully nonlinear and dissipative solitons in media with spatially-inhomogeneous 
nonlinear absorption, but requiring linear gain, have been reported in Il25ll . 


3. Dynamics of ideal Airy beams in the nonlinear medium 

We turn now into the problem of the propagation of Airy beams in nonlinear media. We are 
interested in an intensity regime where the effects of NLLs are significant, being comparable 
to, or even dominant over other nonlinear effects (Kerr nonlinearity in our case). This means 














Fig. 3. For a medium with M = 5 and a = 1, (a) Medial: Change of the intensity profile 
of an input Airy beam with ag = 3.5 with propagation distance v. The red curve represents 
the final attracting NAB with ag = 4.13, whose inward ffankel amplitude is |a; n | = 3.5. (b) 
NLLs N-oo = 2 \A\ 2 du (black curve) and peak intensity (gray curve) of the propagating 
beam in (a) as a function of propagation distance v. (c) For M = 5 and a = 0 (no Kerr 
nonlinearity), in green, amplitude ao(°°) of the attracting NAB (solid curve), and its inward 
and outward ffankel amplitudes, |a in (°o)| and |aout(°°)| (dashed curves), as functions of the 
amplitude Qo(0) of the input Airy beam. The same lines in black correspond to a = 1 (Kerr 
medium). 


in particular that the possible solitary structures emerging from the Airy beam due to the Kerr 
nonlinearity, as described in IfTOl , are significantly affected by NLLs. This condition is fulfilled 
if the intensity ng/(kgri 2 Xg) required for soli ton formation is higher than, or of the order of 
the typical intensity {2 /at which NLLs are significant. In our dimensionless 
variables, this means that a does not exceed unity significantly. 

First, we address the case of an ideal Airy beam A(n,0) = agAi{u) at the entrance plane v = 0 
of the medium. This ideal situation allows us to distinguish the effects that are due to pure 
nonlinearities from those arising from apodization in finite-power Airy beams (which is known 
to cause also loss of contrast 0). From the point of view of numerical simulations of Eq. ©, 
the strongly distorting effects caused by the finiteness of the computational box in unapodized 
Airy beams are eliminated by the procedure of replacing the propagated envelope A at each 
numerical propagation step by the linearly propagated Airy beam in a narrow interval about 
the ends of the computational box. This procedure is justified because we know that nonlinear 
effects in the beam profile start, as for Bessel beams H 1911201 at the Airy main maximum and 
propagate outwards. Of course, the numerical simulations are valid only up to the distance at 
which nonlinear effects reach the end of the computational box. 

Each Airy beam is then numerically seen to undergo a transformation into a particular sta¬ 
tionary NAB. An example of the dynamics toward the final stationary NAB is shown in Fig. 
[3a), where the input linear Airy with ag = 3.5 in a medium with M = 5 and a = 1 transforms 
into the NAB defined by ag = 4.13. Accordingly, the total NLLs experienced during propaga¬ 
tion and the peak intensity, shown in Fig. [3b), are seen to stabilize into a constant value equal 
to the NLLs N-oo and peak intensity of the NAB with ag = 4.13. Figure[3c) summarizes the 
results of extensive numerical simulations. For each value of ag of the input Airy beam at v = 0, 
named ag(0), the value of ag of the attracting NAB, named ao(°°), is represented (solid curves). 
We also depict in Fig. [3 c) the values of |aj n (°°)| and |aout(°°)| of the NAB with ao(°°) (dashed 
curves), as extracted from Fig. [3d)- The crucial observation here is that, irrespective of the 
values of M and a, |aj n (°°)| is given by the identity function |ai n (°°)| = ao(0). Since for the 
input linear Airy beam aj n (0) = a O ut(0) = ag(0), we then infer that the amplitude of the inward 
Hankel component is preserved in the nonlinear propagation. This observation constitutes the 
response to the question of which particular NAB is formed given the input linear Airy beam. 



























The final NAB is that whose inward Hankel amplitude c/ ln equals to the amplitude ao of the 
input Airy beam. 

Here the question of whether a stationary NAB is always reached or not arises. Numerical 
simulations with apodized Airy beams in lf2TI indicate that a stationary state may not be reached 
because of an instability of the final NAB, this instability being observed as “tangential” emis¬ 
sions during propagation. These emissions are incipiently seen in Fig. [3ja), but they do not 
prevent our ideal, unapodized Airy beam from reaching the inward Hankel-preserving NAB. In 
additional simulations with increasing values of of a up to a = 3 (increasing relative strength 
of Kerr nonlinearity), these emissions become sharper. The final stationary NAB continues to 
be formed but much more slowly and therefore at much longer propagation distances. With 
ideal, unapodized Airy beams in media with a > 3, the numerical procedure explained above 
requires such a long propagation distance, and therefore wide transversal computational box, 
that it is difficult to draw any conclusion about whether the Airy beam enters into a unstable 
regime or converges eventually to the NAB. As pointed out in ETft . a stability analysis of NABs 
would help to elucidate this question. In this respect the situation is different with Bessel beams, 
where unstable nonlinear propagation regimes are clearly identified 112111261 . 

4. Dynamics of finite-power Airy beams in the nonlinear medium: abrupt and soft input 

We consider now the nonlinear propagation Airy beams with finite power. As a representative 
case, we take the exponential-Airy beam El 

A(u, v) = a 0 Ai(u + i - 7 v) e ‘^(“> v ) e ^“- v2 / 4 ) e ‘ v ^/ 2 (9) 

generated, for example, by imprinting a cubic phase onto a Gaussian beam and using a 2-/ 
Fourier transform system ge As is well-known, the transversal profile of this beam evolves 
from a nearly Gaussian profile at the far field (v —> ±°°) to the Airy-exponential profile 
cioAi(u) exp(aM) at the “waist” (v = 0). We may then wonder whether its nonlinear propaga¬ 
tion has something to do with that of the ideal, diffraction-free Airy beam, particularly when 
the entrance plane of the medium is placed at large negative v, where the profile approaches a 
Gaussian profile. 

Figure [4| a) shows a typical example of the propagation dynamics of the exponential-Airy 
beam in Eq. <(G]) when the entrance plane of the medium coincides with the waist v = 0. We 
will refer this configuration to as “abrupt” input, since the intensity is maximum at the entrance 
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Fig. 4. For a medium with M = 5 and a = 1 with entrance plane at v = 0 (abrupt input): (a) 
Media2: Change of the intensity profile of an initial exponential-Airy beam with ao = 3.5 
and y = 0.0025 with propagation distance v. The red curve represents the ideal attracting 
NAB with ao = 4.13 (|a; n | = 3.5); (b) NLLs and (c) peak intensity versus propagation 
distance v for y = 0, or ideal case (light gray), y = 0.0025 (dark gray), and y = 0.0075 
(black). 

























Fig. 5. For M = 5 and a = 1 and entrance plane at v = —50 (soft input): (a) Media3: 
Change of the intensity profile of an initial exponential-Airy beam with ay = 3.5 and y = 
0.0025 with propagation distance v. The red curve represents the ideal attracting NAB with 
a o = 4.13 (|oti n | = 3.5); (b) NLLs and (c) peak intensity versus propagation distance v for 
y = 0.0025 (dark gray) and y = 0.0075 (black). The horizontal lines represent the NLLs 
and the peak intensity of the attracting NAB. The dashed gray curves in (c) represent the 
peak intensity of exponential-Airy beams propagating linearly. 


plane and therefore the nonlinear effects are initially strong. A detailed comparison between 
Fig. [4] a) and Fig. [2 a) shows that during the first stages of propagation the apodized Airy beam 
follows the same dynamics of transformation towards the attracting NAB. It cannot be reached, 
however, because the finite power flowing towards positive u is consumed before the slow 
process of NAB formation is concluded. The variation of the NLLs and of the peak intensity 
during propagation [Figs.0tb) and [4] c) ] are seen to reproduce initially those of the ideal case, 
but at lower levels of NLLs and intensities that further decrease down to zero. 

Figure prefers to the propagation of the same exponential-Airy beam as in Fig.[4]except that 
the entrance plane in the medium is at v < 0 sufficiently before the waist for the beam to be 
widespread at low intensity levels (|A| 2 <C 1) and for NLLs to be therefore initially negligible. 
The dynamics observed in Fig. [2 a) is independent of the particular choice of the entrance plane 
v < 0 provided that this ’’soft” input condition is satisfied. The remarkable fact is that, in spite 
of the completely different behavior of the transversal profile, the same law of attraction by 
the stationary NAB that preserves the amplitude of the inward Hankel component presides the 
dynamics. In contrast to abrupt input, the attracting NAB is reached in a broad region about the 
focus, a transversal region whose width depends on the apodization strength. Figures [5jb) and 
HJc) show that the NLLs and peak intensity do not present the ripples of the abrupt case, but 
reach smoothly values close to the NLLs and peak intensity of the attracting NAB. 

Thus, in spite of the variety of situations, the propagation of finite-power Airy beams is 
primarily governed by the transformation towards the same attracting NAB as in the ideal case. 
In turn, this NAB is solely determined by the parameter ao of the input Airy beam, or what 
amounts the same thing, by the maximum intensity ~ 0.287«f, of the input linear Airy beam (or 
the maximum intensity that would reach at the focus for soft input). 

Consider finally the limiting strong nonlinear regime of increasing oq. The attracting NABs 
are then characterized by a {) approaching fl( )max , and, as seen in Sec. 12.21 these NABs tend 
to match the LNAB with ao.max in a transversal region expanding towards u —> — From 
a practical point of view, any exponential-Airy beam smoothly launched in the medium will 
result in the partial formation the LNAB irrespective of the large value of a o, as illustrated in 
the example of Fig. [2 a), since apodization prevents from the observation of any (ever further) 
Hankel tails. 

As seen in Figs. [2c) and [2b) for the peak intensity under soft input conditions, the formed 
NAB or LNAB persists for a propagation distance that increases with ao and weakening 
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Fig. 6. For M = 5 and a = 1 and entrance plane at v = —75 (soft input): (a) intensity 
profile at the focus v = 0 of an input exponential-Airy beam with ao = 7.0 and y = 0.0025 
(black curve), and intensity profile of the LNAB in this medium (red curve), (b) Peak 
intensity versus propagation distance for y = 0.0025 (dark gray) and y= 0.0075 (black). 
The horizontal line represents the peak intensity of the LNAB. The dashed gray curves 
represent the peak intensity of exponential-Airy beams propagating linearly. 
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Fig. 7. Dynamics of a linear Airy beam with ao = 2 entering a nonlinear medium with 
M = 5 and a = 0 at v = 0 and exiting it at v = 10. (a) Media4: Change of the intensity 
profile with propagation distance. The red curve represents the initial lineal Airy, (b) Peak 
intensity versus propagation distance. The dashed vertical lines represent the region filled 
by the nonlinear medium. 


apodization. Observing in these figures that the peak intensity rises initially, and decreases after 
the plateau region of intensity about unity, approximately as that of the corresponding linear 
Airy beam (gray dashed curves) ~ 0.287flQexp(—yv 2 /2), we readily obtain the estimation 

Av ~ 2\j 21n(0.287a5)/y (10) 

(provided that 0.287 «q > 1) for the distance of sustained peak intensity. The length of this 
region can be made arbitrarily large not only by weakening the apodization, as for linear Airy 
beams, but also by increasing the intensity of the input Airy beam. 

5. Conclusions 

We have described the nonlinear propagation of Airy beams in media with significant NLLs due 
to multiphoton absorption. Though we have restricted the expressions and examples to media 
with self-focusing Kerr nonlinearity, similar description is expected to hold for other dispersive 
nonlinearities, or with the inclusion of ionization and plasma effects, as is also the case with 
Bessel beams E). 
























As expected from CD, NABs play an essential role in the nonlinear propagation of Airy 
beams. We have provided with a full asymptotic analysis of NABs. This provides the basis 
for understanding the nonlinear dynamics of the initial Airy beam. This asymptotic analysis 
itself shows that there exist two types of NABs: the NABs described in CD with asymptotic 
linear, unbalanced H link el decay as 1 / (—w) 12 in intensity, and the limiting NAB wich is fully 
nonlinear and infinitely lossy, decaying as slowly as 1/(— u) l ^ 2M ~ 2 \ 

Our numerical simulations with initial ideal or finite-power versions of Airy beams reveals 
that a stationary NAB tends always to be formed, and unveils a conservation law that deter¬ 
mines the particular NAB that attracts the input Airy beam, namely, the NAB that preserves the 
inward Hankel component of the input Airy beam. This conservation law allow to understand 
the seemingly different behaviors that are observed in the dynamics with input ideal or real 
Airy beams. 

This law appears to hold under much more general conditions, and therefore deserves further 
consideration in related problems, as in the nonlinear propagation of Bessel beams. In Figs. 
[TJa) andfTJb), for example, an Airy beam enters and exits a medium with NLLs of finite length. 
The conservation law implies what we can call a dissipative but reversible dynamics: The Airy 
beam is attenuated to a NAB, but once the beam leaves the medium, the initial Airy beam is 
restored, as required by the conservation of the inward Hankel component. It is interesting to 
compare this situation with the counterintuitive situation of the damped oscillations by fric¬ 
tion whose initial amplitude is recovered after friction is switched off. This counterintuitive 
behavior, closely related to the self-healing property of Airy beams (6), is possible because the 
Airy beam can be regarded to be an open system that interchanges continuously energy with an 
energy reservoir. 
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